Interplay of Coulomb and electron-phonon interactions in graphene 
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We consider mutual effect of the electron-phonon and strong Coulomb interactions on each other 
by summing up leading logarithmic corrections via the renormalization group approach. We find 
that the Coulomb interaction enhances electron coupling to the intervalley Ax optical phonons, but 
not to the intravalley E 2 phonons. 



Electron-phonon coupling (EPC) in graphene is 
currently a subject of intense research. Experi- 
mental information on EPC is obtained by Raman 
spectroscopyi^^ and angle-resolved photoemission 
spectroscopy (ARPES). 6 Theoretically, EPC constants 
are usually calculated from density-functional theory 
(DFT),— &2 where the exchange (Fock) term is treated 
in the local density approximation (LDA)ji£ or the gen- 
eralized gradient approximation (GGA)»ii Disagreement 
between the calculated EPC and the ARPES data has 
been pointed out in Ref. [||. 

Also, the ratio between calculated EPC constants for 
different phonon modes disagrees with the ratio of the 
integrated intensities In* /Ig* °f the corresponding two- 
phonon Raman peaks, as noted in Ref. H2I. Namely, in the 
Raman spectrum of graphene two two-phonon peaks are 
seen: the so-called D* peak near 2u>Ax — 2650 cm" 1 , and 
the G* peak near 2u>e 2 = 3250 cm -1 , corresponding to 
scalar A\ phonons from the vicinity of the K point of the 
first Brillouin zone, and to pseudovector E2 phonons from 
the vicinity of the T point, respectively. Experimentally, 
Id*/Ig* ~ 20ji which cannot be reproduced with the 
EPC constants obtained from DFT calculations. 

In this work we consider mutual effect of weak electron- 
phonon and strong Coulomb interactions on each other 
by summing up leading logarithmic corrections via the 
renormalization group (RG) approach, which goes be- 
yond Hartree-Fock approximation. Coulomb interaction 
has been known to be a source of logarithmic renormal- 
izations for Dirac fermionsi 13 i 14 ' 15 Coulomb renormaliza- 
tions in graphene subject to a magnetic field have been 
considered in Ref. Coulomb effect on static disorder 

has been studied in Refs. mum. To the best of our 

knowledge, Coulomb renormalization of EPC has never 
been considered; moreover, at energies higher than the 
phonon frequency, EPC itself is a source of logarithmic 
renormalizations, and has to be included in the RG pro- 
cedure. 

Upon solution of the RG equations we obtain that 
(i) EPC tends to enhance Coulomb interaction, but not 
sufficiently to dominate over the flow to weak coupling, 
found earlier^ (ii) Coulomb interaction enhances the 
EPC only for the scalar A\ phonons, while renormal- 
ization of the coupling to the pseudovector E2 phonons 
is only due to EPC and is relatively weak, in agreement 
with the Raman datai 

We measure the single-electron energies from the Fermi 
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TABLE I: Irreducible representations of the group Ce v and 
their characters. 
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TABLE II: Classification of 4 x 4 hermitian matrices by irre- 
ducible representations of the Ce v group. Matrices joined by 
braces transform through each other under translations. 



level of the undoped (half-filled) graphene. The Fermi 
surface of undoped graphene consists of two points, called 
K and K' . Graphene unit cell contains two atoms, each 
of them has one 7r-orbital, so there are two electronic 
states for each point of the first Brillouin zone (we dis- 
regard the electron spin). Thus, there are exactly four 
electronic states with zero energy. An arbitrary linear 
combination of them is represented by a 4-componcnt 
column vector ifi. States with low energy are obtained by 
including a smooth position dependence tp(r), r = (x, y). 
The low-energy hamiltonian has the Dirac form: 20 



d 2 r^(r) (-wS ■ V)^(r). 



(1) 



We prefer not to give the explicit form of the isospin 
matrices £ = (Y* x ,Y, y ), which depends on the choice of 
basis (specific arrangement of the components in the col- 
umn V0- We only note that all 16 generators of the SU (4) 
group, forming the basis in the space of 4 x 4 hermitian 
matrices, can be classified according to the irreducible 
representations of Cq v , the point group of the graphene 
crystal (Tables Q] and [TTJ) . They can be represented as 
products of two mutually commuting algebras of Pauli 
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FIG. 1: Phonon displacements for E2, A\, and B\ modes. 
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FIG. 2: Analytical expressions of the diagrammatic tech- 
nique. 



matrices E x , £„, E z and A x , A a , A z ) 21 i 22 which fixes their 
algebraic relations. By definition, ExjEj, are the matri- 
ces, diagonal in the K, K' subspace, and transforming 
according to the E\ representation of Cq v . The Fermi 
velocity v w 10 8 cm/s. 

The hamiltonian of the long-range Coulomb interac- 
tion between electrons has the form (hereinafter we imply 
the summation over the spin indices): 



2„/ P( r )p( r ') 



p(r)=^t( r )^( r ), ( 2 ) 



The background dielectric constant of the substrate is 
incorporated into e 2 . 

For low-energy electronic states EPC is efficient if the 
phonon wave vector is close to V, K or K' point. Con- 
sidering only in-plane displacements, we have 4 degrees 
of freedom per unit cell. Consider the T point first. Two 
modes are acoustic, they weakly couple to electrons, and 
are neglected. The other two correspond to £2 (pseu- 
dovector) optical phonons, shown in Fig. [TJ They couple 
to the electronic motion via the iCfC '-diagonal £2 matri- 
ces from Table HIl The K and K' points are related by 
the time reversal symmetry, so the phonon frequencies 
are the same, and one can form real linear combinations 
of the modes from K and K' . They transform accord- 
ing to Ai, B\, A 2 , B2, Ei, E2 representations of Cq v , and 
couple to the electronic motion via corresponding KK'- 
off-diagonal matrices. Linear coupling to A2 and B2 dis- 
placements is forbidden by time-reversal symmetry, and 
coupling to the E\, E2 modes is numerically small. 23 The 
reason for this smallness is that Ei , E 2 displacements do 
not change any C-C bond length; in the tight-binding 
approximation this coupling simply vanishes. Thus, we 
restrict our attention to the E2 modes from the T point 
and Ai and Bi combinations of the modes from the K, K' 
points, shown in Fig. [T] They are the only modes seen in 
the Raman spectra of grapheneJi 2 -^^ 

We take the magnitude of the carbon atom displace- 
ment as the normal coordinate for each mode, denoted 
by Uf,,, fi — x, y, a, b for the four modes, shown in Fig. [T] 
respectively. Upon quantization of the phonon field, 
and the phonon hamiltonian _ff p h are expressed in terms 
of the creation and annihilation operators K u , 6 q u, as 



= ^ b qA1 e^ + h.c. 



(3) 
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FIG. 3: Electron self-energy due to the screened Coulomb 
interaction and the EPC. 



The crystal is assumed to have the area L x L y , and to 
contain N carbon atoms of mass M. The q summation is 
performed as ^ q — * L x L y J d 2 q/(27r) 2 . "h.c." stands for 
hermitian conjugate. The two degenerate E2 modes have 
the frequency loe 2 ~ 0.196 eV. As the Ai and B\ modes 
represent real linear combinations of modes from K, K' 
points, they have the same frequency, u>Ai ~ 0.170 eV, 
and appear with the same coupling constant in the EPC 
hamiltonian: 



H 



EPC 



$(r) 



^F^(r)(AE) A 

A" 

= J d 2 v^{v){F E2 [u x (r)A z T, y — u y (r)A z T, x ] + 

+ F Al {u a (r)A x Z z + u b {r)A y Z z ]} i>(r). (4) 



The coupling constants Fe 2 and Fa ± are not related 
by any symmetry. However, in the tight-binding model 
Fe 2 = Fa x = 3(dto/da), where to is the nearest-neighbor 
coupling matrix element, and a is the bond length. 

As we are interested in energies much higher than the 
temperature T, we set T = 0; still, calculations are much 
more transparent in the Matsubara representation. Elec- 
tron and phonon Green's functions are given by 



G(p,ie) 



{vpf 



D^iuj) = -- 



2u„ 



' 2 +< 



(5) 



their graphical representation is shown in Fig. [2] The 
electronic self-energy is a sum of the Coulomb^ and 
EPCH contributions: E = E ee + £ p '\ shown in Fig. O 



The leading logarithmic term in E ee is given byJ^ 

/did dp 
7T^2 V(<k M G(P - q, « - iw) 
27T (27Tj^ 
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7T arccosg - garccos^ 
/(5) = 1 -7T- + — /„ f{9) = — F==, (7) 
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nAfe 
8v 



■(8) 



A/" = 4 is the number of Dirac species, valley and spin 
degeneracy taken into account. The lower cutoff £ m i„ ~ 
max{up, e}, the upper cutoff ^ ma x ~ v/a is of the order 
of the electronic bandwidth. The logarithmic divergence 
in the Fock self-energy S ee is due to long-distance nature 
of the Coulomb interaction, and thus is not picked up by 
local approximations such as LDA or GGA. The random 
phase approximation (RPA) for V(q, iuj), shown in Fig. [31 
corresponds to expansion of the pre-logarithm coefficient 
to the leading order in the parameter 1/Af = 0.25, as- 
sumed to be small. This is justified better than expansion 
in the dimensionless coupling constant g, obtained with 
the bare coupling 27re 2 /q. Indeed, for Af = 4 we have 
g = (ir/2)(e 2 /v) « 3.4; background dielectric screening 
reduces it to g ~ 1. 

The leading logarithmic asymptotics of TP h is given by 



duj d 2 q F 2 
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(9) 



Here £ mi „ ~ max{e,w M }, £ ma2 , - u/a, and \/27a 2 /4 is 
the area per carbon atom. The dimensionless constants 
A,e 2 , Aai will be treated as small parameters. 

The latter statement deserves some discussion. In prin- 
ciple, one could proceed analogously to the Coulomb 
case: instead of doing the perturbative expansion in A^, 
one could dress the bare phonon propagators by the ap- 
propriate polarization operators n(q, iui), corresponding 
to 1/Af expansion. Since II(q, iuj) oc q at u> <C q <C 1/a.r^ 
the dressed phonon frequency would grow as y/q, and 
would no longer diverge logarithmically. However, 
the inelastic X-ray scattering data for the phonon disper- 
sions* show that the phonon dispersion is smaller than 
the phonon frequency itself. Thus, the renormalization 
of the phonon frequency remains small even at q ~ 1/a, 
so the perturbative expansion in A M is more justified. 

The logarithmically divergent integrals in Eqs. © 
and ([9]) have different structure due to different form of 
the screened interaction V(q, itu) and the phonon propa- 
gator D^iiuj). In Eq. ^ the integral is dominated by the 
frequencies \u\ ~ vq, while in Eq. ^ it is \lu\ ~ u)^, since 
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FIG. 4: Logarithmic corrections to the EPC vertex F M of the 
order 0(l/N,\%). Dia grams (c)-(f) vanish. 

i i 

FIG. 5: Logarithmic diagrams of the order O(A^) not reduced 
to a renormalization of the EPC vertex F M . 

D^(iuj) oc 1/w 2 at \u>\ 3> Thus, in the calculation of 
the leading logarithmic asymptotics it is sufficient to ap- 
proximate D^iiuj) ps — 2ttS(u>). This substitution makes 
the phonon propagator (combined with EPC vertices) 
formally analogous to the correlator of static disorder 
potential (i. e., from the point of view of electrons with 
t Up the lattice is effectively frozen). Thus, renormal- 
izations due to EPC at e > w,, are equivalent to those 
due to static disorde r 17 i 18 i 19 i 22 . This equivalence holds 
only in the leading order in EPC, since in higher orders 
the phonon propagator is dressed by polarization loops, 
and the static disorder correlator is not. 

The presence of the large logarithm invalidates the 
first-order expansion in 1 /Af, and makes it necessary to 
sum all leading logarithmic terms ~ (1/Af) n In™ of the 
perturbation theory. This is done using the standard RG 
procedure i 13 i 14 ' 15 i 16 ' 17 i 18 ' 19 Let us introduce the running 
cutoff £ m a X e~ e - One RG step consists of reducing the 
cutoff, I -► 1+ 81, so that e~ se < 1, while (1/Af)8l < 1, 
X^Si <C 1. The inverse Green's function transforms as 



it — t>p£ — S(p, it) = 



it — (v + (5u)pS 
1 + 8Z 



(10) 



where SZ is chosen to preserve the coefficient at it upon 
rescaling of the electronic fields, t/> — > (1 + 8Z/2)ip: 8Z = 
9S(p, it)/d(it). Then v is renormalized as 



Sv dE(p,it) 9E(p,ie) 



d(it) 



a(wps) 



(ii) 



Next, we determine renormalization of the coupling 
constants. The electron charge e is not changed, as guar- 
anteed by the gauge invariance, so the renormalization 
of the Coulomb coupling constant g is determined by the 
velocity v. For the EPC, logarithmic vertex corrections 
of the order 0(1/ 'Af, A 2 ) are shown in Fig. 0] Two other 
diagrams (Fig. O should be taken into account, as they 
are of the same order and also logarithmic. As a result, 
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FIG. 6: Flow of the dimensionless coupling constants X A i 
(three upper curves, starting from the bare value 0.04 at 
10 eV, respectively) for three values of the bare Coulomb cou- 
pling g(0) — 3.4, 1.5, 0.5 (solid, dashed, and dotted curves). 
The constant Xe 2 = 0.035 is unchanged in the order O(X). 



we obtain the following RG equations: 



dg 
di 
d\E 2 
dl 
d\Ai 
d£ 



8/(g) 



Ag 2 + Aai 
2tt 



16/(ff) 



Ayll ■ 



(12a) 
(12b) 
(12c) 



Because of the diagrams of Fig. El the renormalized A p 
cannot be related to a new EPC vertex F^. Iterations 
of these diagrams generate electron coupling to multi- 
phonon excitations, all included in the renormalized A M . 

The coupling constant Xe 2 at energies ~ u>e 2 ~ 0.2 eV 
can be extracted from the experimental data of Ref. HI 
(change of the Raman G peak with the electron den- 
sity): \e 2 ~ 0.035, corresponding to Fe 2 = 6 eV/A. 
Then the main effect comes from the Coulomb terms: 
as f(g) > 0, g flows to weak coupling,- 1 - Ax a is en- 
hanced, while the enhancement of Xe 2 , proportional to 
X 2 Ai , is much weaker due to the cancellation between 
Coulomb self-energy and vertex corrections. This be- 
havior in qualitative agreement with the Raman data: 



when the A 2 term is neglected, the ratio of the inten- 
sities of the two-phonon peaks, mentioned in the intro- 
duction, is Id'/Ig' = 2(A J 4 1 /A£ 2 ) 2 ji2 If only Coulomb 
terms are kept in Eqs. Q12p , their integration gives 
X A} {t)/X M {Q) = [g(0)/g(£)} 2 = [v(£)/v(0)} 2 , which, in 
principle, can be checked experimentally. 

To study the behavior of the coupling constants quanti- 
tatively, we solve Eqs. (|12p numerically, neglecting the A 2 
term. The largest value of I is determined by the lower 
cutoff ^min ~ ~ 0.2 eV. In Fig. [6j we show the 
flow of \a x for three values of the bare Coulomb cou- 
pling constant: g(0) — 3.4 (corresponding to no dielec- 
tric screening at all), 5(0) = 1.5, and 5(0) = 0.5. The 
bare values of the the electron-phonon coupling constants 
A B2 (0) = 0.035, X Al (0) = 0.040 were chosen (a) to sat- 
isfy the relation Xe 2 (0) / X Al (0) = lu Ai /lue 2 , valid in the 
tight-binding approximation, (b) to reproduce the exper- 
imental value Xe 2 ~ 0.035. Note that the RPA calcula- 
tion without the RG collection of all leading logarith- 
mic terms, would give all dependencies on Fig. [6] to be 
straight lines with slopes fixed at 10 eV. A comparable er- 
ror would be produced by the GW approximation, which 
neglects vertex corrections, and thus picks up correctly 
only the first term of the logarithmic series. 

To estimate the EPC strength relevant for Raman scat- 
tering, we stop the RG flow at electronic energies ~ 1 eV 
(half of the incident laser frequency). In the unscreened 
case, g(0) — 3.4, it gives X Ai /Xe 2 ~ 3.2, in agreement 
with the observed ratio Id* /Ig* ~ 20. 

To conclude, in this paper we have considered mutual 
effect of the weak electron-phonon and strong Coulomb 
interactions on each other by summing up leading log- 
arithmic corrections via the renormalization group ap- 
proach in the intermediate energy range u>e 2 ,u> Ai < 
e < v/a. At these energies quantum fluctuations of the 
phonon field may be viewed as effective static disorder. 
We find that Coulomb interaction enhances electron cou- 
pling to the intervalley A± optical phonons, but not to 
the intavalley E2 phonons, in agreement with the exper- 
imental data for two-phonon Raman scattering. 

We thank M. S. Foster, F. Guinea, and F. Mauri for 
helpful discussions. 
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